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Quantum  Mechanical  Approach  to  Understanding 
Microstructural  and  Mechanical  Properties  of  Intermetallics 


by 

Arthur  J.  Freeman 
PRINCIPAL  INVESTIGATOR 


I.  INTRODUCTION  AND  BACKGROUND  PERSPECTIVE 

The  high  temperature  intermetallics  have  become  a  major  focus  of  efforts  -  both 
experimental  and  theoretical  -  driven  in  large  part  by  their  potential  for  important  aerospace 
applications.  In  the  theoretical  area  of  this  research,  our  woik  has  been  at  the  forefront  and 
has  already  stimulated  similar  efforts  at  other  laboratories.  Indeed,  the  work  initiated  and 
carried  out  under  AFOSR  support  reflects  the  present  state  of  the  field  and  the  reaching  of  a 
threshold:  advanced  theoretical  -  computational  techniques  combined  with  the  power  of 
supercomputers  provide  an  understanding  of  matter  at  the  atomic-scale  with  an 
unprecedented  level  of  detail  and  accuracy.  It  is  this  capability  that  has  given  birth  to  a  new 
branch  of  scientific  endeavor;  Computational  Materials  Science.  Not  only  can  we  simulate 
experiment,  but  using  the  computational  approach,  we  are  now  beginning  to  design  new 
materials  and  to  predict  their  properties  without  actually  synthesizing  them.  In  contrast  to  an 
analytic-theoretical  approach,  which  isolates  and  idealizes  real  systems  to  unravel  fundamental 
relations  and  laws,  the  computational  approach  is  synthetic:  its  goal  is  to  simulate  more  and 
more  details  of  the  systems  studied  including  as  much  of  the  environment  as  possible. 

Designing  new  materials  with  specific  mechanical,  thermal,  chemical,  electronic, 
magnetic  and  optical  properties  hinges  on  one  basic  assumption;  the  properties  of  the 
macroscopic  ensemble  are  related  to  and  can  be  derived  from  the  properties  of  individual 
molecules  and  atomic  building  blocks  (such  as  crystallographic  unit  cells  in  solids).  But 
which  atomic-scale  physical  quantities  and  observables  need  to  be  looked  at?  This  is  the 
important  question  today  and  the  key  to  successful  materials  by  design  in  the  fumre. 

Today  much  of  the  theoretical  efforts  are  focused  on  the  atomic  scale  since  the 
theoretical  concepts  and  methodologies  being  developed  are  bringing  us  closer  to  the  goal  of 
a  unified  approach.  The  theoretical  approach,  density  functional  theory,  that  underlies  all  the 
calculations  discussed  here  is  appealing  because  it  promises  to  provide  such  a  unified 
approach,  which  is  particularly  suited  for  bridging  the  gap  between  our  understanding  of 
atomic  scale  properties  and  the  macroscopic  properties  essential  for  producing  useful 
materials  applications. 
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Our  goal  was  to  study  and  develop  alloying  concepts  for  understanding  intennetallic 
alloys  as  derived  from  a  first  principles  quantum  mechanical  approach.  Thus,  a  major  part  of 
our  effort  was  to  study  and  determine  effects  of  alloying  elements  in  several  intermetallic 
compounds  and  to  work  closely  with  experimental  efforts  to  evaluate  the  applicability  of  the 
theoretical  approach  to  alloy  design.  Specifically,  as  seen  from  the  publications  listed,  our 
highly  precise  all-electron  quantum  mechanical  electronic  stracture  methods  were  applied  to 
the  study  of  a  number  of  materials  problems  in  order  to  obtain  from  first  principles, 
information  of  relevance  to  alloy  stability  and  the  design  of  structural  materials.  Using  our 
recently  developed  state-of-the-art  all-electron  self-consistent  total  energy  methods  which 
give  precise  solutions  of  the  local  density  equations,  fundamental  information  was  obtained 
about  the  structural  and  electronic  properties  of  these  alloys  in  order  to  predict  stable  and 
metastable  phases  and  how  alloying  affects  bonding,  crystal  ordering  and  crystal  symmetry. 
The  first  principles  local  density  approach  used  addresses  questions  of  a  metallurgical  nature, 
such  as  phase  stability,  crystal  structure,  equilibrium  lattice  constants,  and  some  mechanical 
properties  including  the  effect  of  atomic  relaxation. 

We  have  attempted  in  the  above  to  describe  how  our  first  principles  approaches  have 
been  employed  to  study  some  key  problems  in  understanding  high  temperature  intermetallics. 
It  needs  to  be  understood  that  the  problems  posed  are  as  con^lex  and  challenging  as  any  in 
materials  science  today.  Thus,  while  significant,  the  progress  made  by  us  and  others 
constitute  (of  necessity)  only  the  first  steps  taken  to  "bridge  the  gap"  -  see  the  accompanying 
table  (developed  by  Rosenstein  &  Freeman)  -  of  relating  the  macroscopic  properties  to  the 
atomic  scale  properties.  This  table  encompasses  not  only  a  philosophy  but  also  a  working 
principle,  that  bridging  the  gap  requires  a  strong  collaboration  between  the  materials  science 
and  metallurgy  community  (Dimiduk,  Frazer  and  others)  and  computational  materials 
scientists. 

Since  most  of  our  recent  effort  focussed  on  y-TiAl,  we  describe  some  of  this  work 

here.  This  will  also  serve  to  illustrate  the  approach  taken  in  tackling  some  complex  problems, 
the  computational/theoretical  capabilities  brought  to  bear  in  their  solution  and  the  meaningful 
results  obtained. 

n.  PROGRESS  MADE 

A.  Electronic  Structure,  Ti-Ti  Bonding  and  Properties  of  y-TiAl 

Recently,  several  authors^2,l3,14  have  tried  to  understand  the  cause  of  the  poor 
ductility  of  y-TiAl  in  terms  of  the  type  and  character  of  available  deformation  mechanisms. 

Some^5,l6,6,l7  have  tried  to  examine  possible  effects  of  bonding  on  mechanical  properties 


3 


COMPUTATIONAL  MATERIALS  DESIGN:  HIGH  TEMPERATURE 
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and  to  understand  the  dislocation  behavior  in  terms  of  the  underlying  bulk  electronic 
structure.  Studying  the  effects  of  substitutional  ternary  additions  using  a  cluster  linear 

muffin-tin  orbital  (LMTO)  method,  Morinaga  et  al.l5  observed  that  ternary  additions  (such  as 
Mn)  which  weaken  tiie  Ti-Al  d-p  interactions  and  enhance  the  H-H  d-d  interactions  may 
improve  its  ductility.  Using  the  full-potential  linearized  augmented  plane  wave  (FLAPW) 

method^  Fu  and  Yoo  studied  the  elastic  constants  and  various  fault  energies  of  TiAl  and 
found  a  fairly  large  elastic  shear  anisotropy  (A2)  along  the  [Oil]  direction,  and  a  high  APB 
energy  due  to  the  strong  cohesion  between  'll  and  A1  layers.  From  charge  densities  and 
planar  fault  energies  calculated  using  a  layered  Korringa-Kohn-Rostoker  method.  Woodward 

et  al.^^  observed  that  the  hierarchy  of  fault  energies  was  consistent  with  the  number  of 

disrupted  p-d  and  d-d  bonds.  Greenberg  et  al.^  suggested  that  'll  d-d  bonding  in  the  (001) 
plane  forms  a  rigid  framework  and  that  this  charge  anisotropy  results  in  a  Peierls  stress 

anisotropy  in  y-TiAl.  In  their  view,  the  presence  of  deep  Peierls  valleys  may  be  one  factor 

responsible  for  both  brittleness  and  the  anomalous  flow-stress  versus  temperature  behavior. 

It  is  well-known  that  the  flow  and  fracture  of  high  temperature  intermetallic  alloys  is 
strongly  influenced  by  the  propagation  and  cross  slip  of  ordinary  and  super  dislocations.  In 
TiAl,  there  is  an  apparently  large  Peierls  stress,  or  lattice  friction  stress  for  the  glide  of 
ordinary  (1/2<110>)  dislocations.  This  Peierls  relief  is  intimately  related  to  the  underlying 
crystal  and  electronic  structure  of  TiAl,  and  specific  bonding  states  have  been  identified  by 
some  investigators  as  contributing  to  the  pinning  of  dislocations  along  certain  line  directions. 
In  our  study,  we  quantify  the  relative  strength  of  these  bonds  using  several  complementary 
electronic  structure  methods. 

Many  factors  have  been  found  to  affect  the  ductility  of  this  material,  such  as  ternary 
additions,  microstructure,  and  especially  stoichiometry.  For  example,  the  ductility  of  y-TLAl 

increases  with  (i)  decreasing  aluminum  concentration  (i.e.,  depends  on  the  stoichiometry), 
(ii)  decreasing  amounts  of  impurity  concentration  (such  as,  oxygen  or  carbon  etc.),  (iii) 
doping  few  percent  ternary  additions  ( such  as,  Mn,  Cr,  and  V  etc.).  Many  speculations  were 
proposed  to  discover  a  mechanism  for  improving  its  mechanical  properties.  One  such 
conjecture  is  that  the  increase  in  metallic  bonding  with  decreasing  aluminum  content  may  be 
responsible  for  improving  the  ductility. 

In  what  follows,  aspects  of  the  electronic  structure,  equilibrium  lattice  constants  and 
bond  energies  of  y-TiAl  are  compared  with  hypothetical  fee  U  as  obtained  by  the  full- 

potential  linearized  augmented  plane  wave  (FLAPW)  method^,  the  DMol  molecular  cluster 
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method^  and  the  linear  muffin-tin  orbital-Green’s  function  (LMTO-GF)  method^ .  Pair 
energies,  bond  occupations  of  Ti-H  and  Ti-Al  bonds  in  Llo  TiAl  and  fee  H  are  calculated 
using  the  LMTO-GF  method.  Comparing  the  dxy  bonding  states  for  these  two  crystal 
structures,  we  find  that  the  in-plane  directional  d-d  bonding  between  Ti  atoms  on  the  (001) 
plane  are  strengthened  in  TiAl  relative  to  fee  Ti.  These  observations  suggest  that  the 
alternating  (001)  planes  of  A1  and  Ti  atoms  enhance  these  bonding  states.  Conversely, 
transition  metal  ternary  additions  substituted  on  the  A1  sites  in  TiAl  arc  expected  to  weaken 
these  bonding  states.  Using  the  FLAPW  method,  we  examined  the  changes  in  charge  density 
in  the-  (001)  Ti  planes  when  Mn  atoms  are  substituted  on  A1  sites  in  TiAl  and  LMTO-GF 
methods  were  used  to  study  the  effect  of  Mn  substitutional  point  defects  on  the  (001)  Ti  dxy 
bonding  states.  The  elastic  constants  of  TiAl  and  A1  rich  TiAl  were  determined  using  the 
FLMTO  method  and  compared  with  those  of  fee  Ti.  Most  recently,  the  effects  of  ternary 

additions  on  the  twin  energy  and  site  preference  in  y-TiAl  were  investigated  as  was  the  role  of 

impurities  such  as  C  and  O  (see  below). 

B.  Effect  of  Stoichiometry 

It  is  known^2  j^at  y-TiAl  has  the  ordered  face-centered  tetragonal  LIq  structure  with 

c/a=  1 .02,  and  that  single  phase  y-TiAl  exists  only  on  the  Al-rich  side  with  a  composition 

range  about  50  -  65  at.%Al.  At  off-stoichiometric  compositions,  excess  titanium  and 
aluminum  atoms  occupy  anti-sites  without  creating  vacancies.  On  the  other  hand,  the  phase 
diagram  shows  that  Ti-rich  TiAl  decomposes  into  a  mixture  of  two  phases,  TiAl  and  TisAl 

(a2/y). 

For  studying  the  effect  of  stoichiometry  theoretically,  we  constructed  a  16-atom  LIq 
stmetured  supercell  as  a  model  to  simulate  the  Al-rich  (or  Ti-rich)  TiAl  according  to  the  center 
site  occupied  by  A1  (or  Ti).  In  order  to  judge  the  phase  stability,  we  calculated  the  formation 

energy^^  for  each  phase,  and  compared  with  that  of  its  decomposed  components  or  existing 
equilibrium  compounds  (e.g.,  we  compared  the  formation  energy  of  Ti7Al9  with  that  of 
7TiAl-i-2Al  or  6TiAl+TiAl3).  The  calculated  formation  energy  results  show  that  Al-rich  TiAl 
(i.e.,  Ti7Al9)  is  energetically  stable  as  compared  to  its  decomposed  counterparts  (either 
7TiAI-i-2Al  or  6TiAl-i-TiAl3).  On  the  other  hand,  Ti9Al7  is  unstable  with  respect  to  a  mixture 
phases  of  TiAl  and  T13AI  (cf.,  Table.  1  of  our  as  yet  unpublished  results)  in  agreement  with 

observations. 
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Table  I.  Equilibrium  lattice  constant  (in  A)  and  formation  energy  (in  eV/atom)  of  Al-rich  TiAl 
and  Ti-rich  TiAl.  As  comparison,  the  formation  energy  of  the  decomposed  components  or 
existing  stable  compounds  were  also  listed  (unpublished). 


T  attify.  Cnnstant  Formation  Energy 


Ti7Al9 

3.97 

0.649 

7TiAl-h2Al 

0.472 

6TiAl+TiAl3 

0.509 

Ti9Al7 

3.98 

0.340 

6TiAl-hTi3Al 

0.477 

C.  Electronic  Structure  of  TiAl 

In  order  to  gain  insight  at  the  microscopic  level  into  the  bonding  character  of  TiAl,  we 
considered  tetragonal  LIq  TiAl  along  with  hypothetical  fee  Ti  because,  except  for  c/a=1.02, 
tetragonal  face-centered  LIq  is  indeed  a  fcc-like  structure.  One  can  view  LIq  TiAl  as  stacking 

alternatively  one  fee  H  layer  on  one  fee  A1  layer.  To  understand  the  bonding  character  of  y- 
TiAl,  we  studied  the  electronic  structures  of  both  LIq  TiAl  and  fee  Ti  using  the  FLAPW 
method,  and  obtained  the  equilibrium  lattice  constants  for  both  fee  Ti  and  LIq  TiAl.  The 

calculated  lattice  constant  for  y-TiAl  agrees  well  with  the  observed  value.  It  is  interesting  to 

note  that  the  calculated  interatomic  distance  (i.e.,  bond  length)  of  the  Tl-'Ii  bond  in  the  fee 
structure  is  quite  close  to  the  observed  bond  length  of  H  in  the  bcc  structure,  and  about  3% 
smaller  than  that  of  Ti  in  the  hep  structure.  However,  it  is  1%  larger  than  the  li-Ti  bond 
length  in  TiAl.  This  means  that  when  one  'll  ((X)l)  layer  and  one  A1((X)1)  layer  are 
successively  stacked  together  along  the  <(X)1>  direction  to  form  TiAl,  it  keeps  the  fcc-like 
structure;  the  only  modification  seen  is  the  appearance  of  a  shrinkage  of  its  lattice  constant  (or 
the  bond  length)  in  the  ((X)l)  plane. 

To  visualize  the  bonding  character,  we  plotted  the  charge  density  of  Ti-plane  in  the 
pure  fee  Ti,  compared  it  with  its  counterpart  in  LIq  'HAl  and  found  that  the  charge  becomes 
more  directional  (dxy)  in  TiAl  as  compared  with  that  in  fee  Ti.  This  implies  that  the  existence 
of  A1  in  the  alternating  plane  in  TiAl  may  be  responsible  for  this  more  directional  bonding. 
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D.  Bond  Energy  and  Occupation  of  the  Ti-Ti  bonding  in  TiAi 

In  order  to  have  a  quantitative  description  of  the  bonding  character,  we  calculated  pair 
energies  and  occupation  numbers  of  the  bond  with  the  LMTO-Green  Function  method.  The 
in-plane  Ti-Ti  pair  energy  increases  from  fee  Ti,  to  LIq  Ti,  to  LIq  TiAl.  This  means  that  the 
bond  strength  between  the  nearest  neighbor  Ti  atoms  increases  in  going  from  fee  Ti  to  TiAl. 
On  the  other  hand,  the  inter-plane  nearest  neighbor  bond  strength  remains  approximately  the 
same  value  for  all  three  cases.  This  is  consistent  with  the  shrinkage  of  the  in-plane  lattice 
constant  for  TiAl  compared  with  fee  Ti. 

Furthermore,  by  studying  the  bond  energy  and  occupation  for  different  orbitals,  we 
found  that  in  the  Ti-Ti  in-plane  bond,  the  bonding  state  with  the  largest  number  of  electrons 
is  the  dxy-dxy  bond.  The  occupation  number  of  the  dxy-dxy  bond  increases  from  fee  Ti,  to 
Llo  Ti  to  in  TiAl,  and  the  bond  energy  increases  from  fee  Ti,  to  LIq  Ti,  to  LIq  TiAl.  It  is 
obvious  that  there  are  more  electrons  and  a  higher  bond  energy  in  the  Tidxy-Tidxy  bond  in 
TiAl  compared  with  fee  Ti.  These  first  bond  strengths  reported  for  TiAl  provide  strong 
evidence  to  support  the  argument  that  there  is  a  strengthened  Tidxy-Tidxy  bond  in  TLA.1. 

E.  Effect  of  Mn  Addition  on  Ti-Ti  Bonding  in  TiAl 

We  also  considered  the  effect  of  ternary  additions  on  the  in-plane  d-d  bonding  by 
replacing  one  of  the  atoms  in  the  A1  layer  by  a  Mn  atom.  The  charge  density  shows  that  fewer 
electrons  are  accumulated  in  the  Ti-Ti  bond  and  that  the  in-plane  d-d  bonding  is  less 
directional.  With  the  LMTO-GF  calculations,  the  Ti-Ti  pair  energy  in  TigAlyMn  is  found  to 
be  smaller  than  that  in  TiAl;  the  Ti-Al  pair  energy  in  TigAlvMn  (Ti-43.75  at  %  A1  -  6.25% 
Mn)  remains  close  to  the  Ti-Al  pair  energy  in  TiAl;  the  H  dxy-Ti  dxy  bonding  state  in 
TigAhMn  has  fewer  electrons  and  a  lower  energy  than  in  TiAl.  These  results  show  that  the 
in -plane  d-d  bonding  is  weaker  and  less  directional  when  Mn  substimtes  A1  in  TiAl. 

Noting  the  experimental  results  that  Mn  additions  occupy  A1  sites^^  and  improves  the 
ductility^O  of  Ti  rich  y-TiAl,  our  results  are  consistent  with  previous  models  that  suggest  that 
the  strengthened  in-plane  directional  bonding  in  TiAl  is  among  the  factors  responsible  for  the 
poor  ductility  of  y-TiAl. 

This  finding  leads  us  to  expect  that  any  ternary  additions  substituted  on  the  A1  sites 
and  contributing  d-electrons  may  reduce  the  in-plane  Ti-Ti  directional  bonding  and  reduce  the 
Peierls  anisotropy  in  TiAl.  Indeed,  this  view  may  explain  the  beneficial  effect  of  some  other 

d-additions  (such  as  Cr  and  V)  in  y-TiAl. 
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F.  Phase  Stability,  Bonding  Character  and  Elastic  Constants  of  TiAl  and  A1  rich  TiAl 
For  a  structural  material,  it  is  important  to  study  its  elastic  constants,  because  the 
elastic  constants  give  a  direct  measure  of  the  stiffness  of  the  materials.  As  was  done  for  the 
determination  of  the  equilibrium  lattice  constants,  the  elastic  constants  can  be  obtained  by 
means  of  total  energy  calculations  under  several  uniform  distortions.  Using  the  FLMTO 

method^,  we  determined  and  compared  the  elastic  constants  for  fee  Ti,  TiAl  and  Ti7Al9. 
These  unpublished  results  are  listed  in  Table  2,  which  is  taken  from  Appendix  A,  the  preprint 
of  a  paper  scheduled  to  appear  in  Phys.  Rev.  B. 

As  shown  in  Table  2,  except  for  C44.  C66  and  C66,  all  the  other  elastic  constants  for 
TiAl  and  the  Al-rich  compound,  Ti7Al9,  remain  more  or  less  the  same  as  that  of  pure  fee  Ti. 
On  the  other  hand,  C44,  which  corresponds  to  a  [010]  shear  strain  on  the  (001)  plane,  has 

been  dramatically  enhanced  -  from  0.73  x  10^^  in  fee  Ti,  to  1.21  x  lO^^  N/M^  in  TiAl, 
to  1.12  X  lOll  N/M2  in  Ti7Al9  .  Note  that  0*66  is  the  only  one  that  increases  monotonically 
upon  going  from  fee  Ti,  to  TiAl  to  Ti7Al9.  Since  C'66=(Cii-Ci2)/2  corresponds  to  the  [110] 
shear  strain  on  the  (1 10)  plane,  this  indicates  that  intrinsic  hardness  increases  from  fee  Ti,  to 
TiAl  toTi7Al9,  which  may  possibly  be  correlated  with  the  lack  of  <110]  dislocations  in  Al- 

rich  Ti48Al52.^^ 


Table  2.  Elastic  constants  (in  unit  of  10"  N/M^)  of  fee  Ti,  TiAl  and  T^Al,  calculated  with  the 
FLMTO  method.  The  percent  number  in  the  parentheses  are  the  difference  from  the  TlAl 
results  listed  in  the  second  column,  the  last  column  lists  the  FLAPW  results  by  Fu  et  al.  for 
comparison. 


Elastic  Constants 

fee  Ti 

TiAl 

Ti^Al, 

TiAl  (FLAPW,  Ref.  15) 

c„ 

1.67(-5%) 

1.76 

1.97(+12%) 

1.90(+8%) 

Cn 

1.01(-h2%) 

0.99 

0.93(-6%) 

1.05(-h6%) 

C,3 

1.01  (+3%) 

0.98 

0.89(-9%) 

0.90(-8%) 

C33 

1.67(-7%) 

1.79 

1.75(-2%) 

1.85(+3%) 

C44 

0.73(-40%) 

1.21 

1.12(-7%) 

1.20(-1%) 

C66 

0.73(+22%) 

0.60 

0.87(+45%) 

0.50(-17%) 

C’s, 

0.33(+15%) 

0.39 

0.52(-h33%) 

0.42(+8%) 
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G.  Effects  of  ternary  additions  on  the  twin  energy  in  y-TiAl 

It  has  been  known  that  in  the  Ti-48%  A1  alloys  the  additions  of  Mn,  V  and  Cr  results 

in  increased  twinning  activity,  and  improves  the  ductility^O.  We  have  attempted  to  understand 
this  behavior  from  a  microscopic  point  of  view,  and  studied  the  (111)  twin  boundary  energy 
of  stoichiometric  TiAl,  of  TiAl  with  ternary  additions  (Mn  and  V),  and  the  effect  of 

stoichiometry  on  the  twin  energy22.  calculations  were  performed  using  a  supercell 

scheme  and  the  total  energy  linear  muffin-tin  orbital  method.  Both  Al-iich  and  Ti-rich  y-TlAl 

were  simulated  with  the  use  of  Ti5Al7  and  TiyAls  supercells,  respectively.  While  the 
calculated  results  presented  may  be  viewed  as  only  qualitative,  the  trends  observed  are 
probably  correct. 

The  calculated  results  showed  that  the  twin  energy  for  y-TiAl  is  about  9  -  10 

mRy/cell,  (equivalently,  70  -  80  erg/cm^),  which  is  comparable  with  the  FLAPW  result^^  (60 
erg/cnP').  Note  that  from  the  twin  energy,  y-TiAl  is  close  to  the  category  of  materials  such  as, 

aluminium,  that  have  high  twin  energy  (75  erg/cm2)23.  The  fact  that  only  few  twins  are 
observed  in  Al-rich  y-TLAl  may  be  understood  from  its  relatively  high  twin  energy.  Indeed, 

for  the  Al-rich  case  the  nearly  same  value  (75  -  80  crg/cnfi)  of  the  twin  energy  is  obtained  as 
that  for  stoichiometric  y-TLAl.  On  the  other  hand,  a  significantly  reduced  twin  energy  (60 

erg/cm2)  was  found  for  the  Ti-rich  case  as  compared  to  stoichiometric  y-TiAl.  This  appears  to 
agree  qualitatively  with  the  observation  of  increasing  twin  activity  when  the  A1  content  is 

decreased  to  48%  in  TiAl^^,  but  note  that  we  have  not  considered  the  mobility  of  twinning 
dislocations  which  should  be  influenced  by  solutes. 

For  studying  the  effects  of  ternary  additions  on  the  twin  energy  in  y-TiAl,  we 
constructed  (Ti5X)Al6  (Tiso-xXxAlso)  and  T16AI5X  (TisoAlso-x)  (where  X  =  Mn  or  V) 
supercells,  representing  V  substitutions  on  Ti  and  A1  sites.22  For  the  case  of  substitutions  on 

the  A1  sites,  the  calculated  twin  energies  turn  out  to  be  55  -  60  erg/cm^,  which  is  about  the 
same  value  as  the  twin  energy  of  Ti-rich  TiAl  but  markedly  lower  than  that  of  stoichiometric 

y-TiAl.  The  observation  of  the  increase  of  twin  activity  in  Ti-48%  A1  manganese  alloys^^ 
may  be  understood  in  terms  of  lowering  the  twin  energy  due  to  the  presence  of  these  ternary 
additions  in  y-TiAl.  On  the  other  hand,  when  V  is  substituted  on  the  Ti  sites  the  calculated 
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twin  energy  is  slightly  lower  than  that  of  pure  y-TiAl.  When  Mn  occupies  H  sites,  an 

exceptionally  small  twin  energy  is  obtained. 

In  order  to  understand  this  from  the  microscopic  point  of  view,  we  inspected  the 
density  of  states  (DOS)  for  Ti6Al5Mn  (or  V).  Near  the  interface,  the  characteristic  p-d  peak 

located  at  about  2.0  eV  below  Ep  in  stoichiometric  y-'KAll^  was  found  to  be  spread  out 

which  indicates  that  the  degree  of  p-d  (directional)  bonding  has  been  reduced.  As  a  result,  the 
general  feature  of  the  total  DOS  for  Ti^AlsMn  (or  V)  resembles  pure  fee  titanium  (in  the 
region  between  Ep  and  -3  eV)  plus  a  parabolic-like  tail  (from  -3  eV  to  -9  eV  (mainly)  due  to 
A1  contributions).  From  the  calculated  partial  DOS,  it  is  clear  that  the  main  contributions  from 
the  ternary  addition  are  the  d-states  located  near  (or  above)  Ep,  which  are  hybridized  strongly 
with  the  Ti-d  states.  This  hybridization  determines  the  major  feature  of  the  DOS  near  Ep,  and 
restores  the  d-d  interactions  between  the  Ti  layer  and  flie  neighboring  d-substituted  A1  layer. 
As  a  result,  the  in-plane  Ti  d-H  d  directional  bonding  is  weakened,  as  found  also  above. 
Thus,  both  reduced  p-d  bonding  and  restored  d-d  bonding  between  neighboring  layers  may 
play  a  positive  role  towards  improving  the  intrinsic  ductility.  (We  should  note  that  a  similar 

analysis  for  the  ylal  and  y/y  interfaces  proposed  later,  will  yield  important  information  there 
as  well.) 

H.  Molecular  Cluster  Studies  of  the  Effect  of  Impurities  in  TiAl 

A  number  of  recent  experiments  show  that  the  mechanical  properties  of  y-TiAl  and 

TiAVTisAl  lamellar  vary  under  different  environments  such  as  in  O2,  H2,  air  or  vacuum.  The 
cause  of  these  environmental  effects  is  still,  however,  unclear.  The  molecular  cluster  model 
represents  a  convenient  method  for  studying  those  electronic  properties  that  are  primarily  a 
function  of  the  local  environment  of  the  system,  such  as  vacancies,  interstitial  and 
substitutional  defects,  and  local  moments.  Thus,  we  have  carried  out  DMol  molecular  cluster 

calculations^,  to  investigate  the  role  of  these  low  level  impurities  on  the  properties  of  TiAl. 
With  this  approach,  the  site  preference,  lattice  relaxation,  and  bonding  character  can  be 
obtained  for  these  interstitial  impurities. 

From  these  calculations,  results  such  as  total  binding  energy  of  the  cluster,  the  charge 
density,  the  deformation  density,  molecular  orbitals  and  density  of  states  were  obtained. 
When  determined  for  several  geometries  or  for  different  electronic  states,  this  information  can 
be  used  to  predict  such  properties  as  relative  energetics  and  stabilities,  and  ground  state 
electron  configurations.  With  atomic  force  calculations,  we  also  determined  the  impurity- 
induced  lattice  relaxation,  which  is  an  important  part  of  the  study. 
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We  began  by  studying  the  electronic  structure  of  TiAl  using  the  total  energy  and 
atomic  force  capabilities  of  the  DMol  method  with  several  cluster  models  ranging  from  6  to 
47  atoms.  We  found  that  clusters  as  small  as  6  atoms  can  give  the  basic  properties  such  as  the 
bonding  between  atoms  -  even  though  there  exists  a  significant  surface  effect.  More  precise 
calculations  should  include  atoms  at  least  to  the  second  nearest  neighbors.  The  results  show 
that  in  the  (001)  pure  Ti  plane,  the  Ti-Ti  directional  bonding  dominates,  which  is  mainly  the 
contribution  of  the  Ti-dxy  component.  In  the  pure  (001)  A1  plane,  the  charge  density  around 
A1  sites  is  quite  flat,  indicating  metallic  bonding  between  the  A1  atoms;  in  the  (010)  mixed  Ti¬ 
Al  plane,  the  Ti-d  and  Al-pz  forms  a  mixture  of  s  and  p  bonding  states.  This  p-d  bonding 
state  gives  a  high  peak  in  the  DOS  curve  at  2  eV  below  Ep.  These  results  are  in  good 

agreement  with  other  theoretical  calculations  and  those  reported  above. 

To  obtain  insight  into  how  the  Ti-Ti  bonding  is  changed  from  the  pure  Ti  metal  to  the 
TLAI  compound,  we  constructed  a  Ti4o  cluster  with  fee  stracture  and  calculated  the 
deformation  charge  between  Ti  atoms  in  the  (001)  plane.  It  turns  out  that  there  is  no  obvious 
directional  bonding  between  Ti  atoms  in  fee  Ti,  in  contrast  with  that  in  TiAl.  This  again 
implies  that  A1  may  play  a  role  in  enhancing  the  directionality  of  Ti-Ti  bonding  in  TiAl. 
Accordingly,  one  can  expect  that  adding  more  A1  to  TiAl  may  further  enhance  the  directional 
in-plane  Ti-Ti  bonding.  Our  calculations  on  an  Al-rich  cluster,  do  show  that  Ti-Ti  in-plane 
bonding  is  stronger  than  that  in  a  Ti-rich  cluster.  The  increase  in  Ti-Ti  in-plane  bonding 
results  in  a  shrinkage  of  the  lattice  constant  in  the  basal  plane  -  which  explains  the 
observation  that  the  lattice  constant  of  the  basal  plane  in  TlAl  is  smaller  than  that  of  pure  fee 
Ti  metal.  This  was  separately  confirmed  by  our  FLAPW  total  energy  calculations  (reported 
above)  which  showed  that  the  horizontal  lattice  constant  of  TiAl  shrinks  by  1%  compared  to 

fee  Ti.  These  results  support  the  idea  of  Greenberg  et  al.^  that  Ti-Ti  bonding  is  also  an 
important  factor  in  determining  the  stmeture  of  TiAl. 

Knowing  the  basic  electronic  structure  of  pure  TiAl,  we  next  began  to  study  the  effect 
of  impurities  on  the  atomic  and  electronic  structure  of  TiAl  and  to  determine  the  possible 
connection  between  these  effects  and  the  mechanical  properties  of  TiAl.  We  studied  two 
impurities,  namely,  carbon  and  oxygen  using  two  clusters  models,  Ti4Al2Xi  and 
Ti2oAl2oXi,  with  X  representing  the  impurity  atom  at  the  octahedral  center  of  TiAl.  There  are 
two  possible  choices  for  the  impurity  atom  position.  One  is  in  the  Ti  plane,  another  in  the  A1 
plane.  Our  calculation  on  the  40  atom  cluster  shows  that  both  C  and  O  prefer  to  stay  in  the  'll 
plane,  since  the  binding  energy  of  this  site  is  more  than  1.9  eV  higher  than  when  C  and  O 
stay  in  the  A1  plane.  Table  3  lists  the  calculated  (as  yet  unpublished)  results.  It  can  be  seen 
from  this  Table  that  C  has  a  higher  binding  energy  than  O  indicating  a  stronger  interaction 
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between  C  and  its  neighboring  atoms.  Without  considering  the  lattice  relaxation,  C  and  O 
have  very  little  hybridization  with  either  Ti  and  A1  implying  that  they  do  not  form  visible 
covalent  bonding  with  metal  atoms.  This  can  be  seen  from  both  density  of  states  and 
deformation  charge  contour  plots,  separately  obtained. 


Table  3:  Binding  energies  (in  eV)  of  40  atom  cluster  with  and  without  impurity  atom. 
XiTi2oAl20  denotes  an  impurity  in  the  Ti  plane;  and  XiAl2oTi20  denotes  an  impurity  in  the  A1 
plane. 


XiTi2oAl20 

Ti2oAl20 

AE 

x=c 

186.91 

177.24 

9.67 

x=o 

185.50 

177.24 

8.26 

XiAl2oTi20 

Al20Ti20 

AE 

x=c 

184.91 

177.40 

7.51 

x=o 

183.75 

177.40 

6.35 

Table  4:  Relaxation  of  the  nearest  neighbor  atoms  when  impurities  are  placed  at  octahedral 

sites  in  the  Ti  plane.  E  (in  eV)  is  the  binding  energy  and  AE  (in  eV)  is  the  energy  gain  after 
relaxation.  Px-Ti  and  Px-Al  stand  for  the  percentage  change  of  X-Ti  and  X-Al  bond  length. 


XiTi20Al20 

E(relaxed) 

E(unrelaxed) 

AE 

Px-Ti 

PX-Al 

X=C 

187.22 

186.91 

0.31 

3.5% 

x=o 

185.96 

0.46 

3.5% 

4.2% 
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An  important  result  from  the  calculation  is  that  we  found  a  significant  relaxation  for 
both  X-Ti  and  X-Al  bond  lengths  when  O  is  placed  on  the  Ti  plane,  while  C  only  induces 
relaxation  for  the  X-Ti  bond  length.  These  results  are  listed  in  Table  4. 

These  results  show  that  O  produces  large  relaxations  in  both  the  H  and  A1  bond 
length  whereas  C  interacts  much  more  strongly  with  Ti  than  with  A1  atoms.  From  the 
difference  charge  density  in  the  X-Ti  plane,  it  can  be  seen  that  there  is  an  obvious  bonding 
between  C  and  the  nearest  Ti  atoms.  For  O  in  TiAl,  this  bonding  is  weak. 

ni.  MOST  RECENT  ACTIVITIES  ON  FIRST  PRINCIPLES  REAL-SPACE  CALCU¬ 
LATIONS  OF  FAULT  ENERGETICS  (INCLUDING  IMPURITIES)  IN  TiAI 

In  this  part  of  the  work,  we  are  starting  to  develop  real-space  methods  for  calculating 
fault  energies.  As  spelled  out  below,  this  is  being  attempted  by  combining  the  first-principles 

tight-binding  Hamiltonian24,25  ^vith  an  effective  real-space  method  for  calculating  the  local 
D0S26.  This  will  hopefully  provide  the  ability  to  study  systems  with  strongly  reduced 

symmetry,  even  amorphous  metals  and  alloys27.28  and  also  impurity  states^^. 

Accurate  ab-initio  calculations  of  stacking  fault  energetics  are  difficult  from  a 

computational  point  of  view.  In  supercell  calculations^®*^  the  numerical  work  increases 

with  the  third  power  of  the  numbers  of  layers;  in  layer  KKR^^.l?  the  planar  defects  are 
ffeated  in  a  less  costly  way;  however,  such  calculations  are  restricted  only  to  planar  defects. 
Moreover,  the  numerical  accuracy  of  these  methods  is  challenged  by  attempting  to  calculate 
small  fault  energies  by  taking  the  difference  of  two  large  total  energies. 

For  systematic  calculations  we  need  methods  which  simplify  the  computational 
difficulties  but  that  are  stiU  ab-initio  and  precise.  Thus,  we  need  to  find  a  compromise 
between  our  wishes  to  simulate  complicated  defects  (such  as  dislocations)  and  to  carry  out 
precise  ab-initio  band  stracture  calculations.  From  our  investigations  to  date,  the  combination 

of  the  ab-initio  tight-binding  like  method  (TB-LMTO)24.25  with  the  recursion  method26  (we 
use  the  notation  R-TBLMTO  method)  and  the  "frozen  potential"  approximation23.34  or  the 

Harris-Foulks  functionaP^  appears  to  provide  a  reliable  approach  for  the  systematic  study  of 
systems  with  reduced  symmetry,  such  as  stacking  faults  in  doped  aluminides.  The 
advantages  of  the  TB-LMTO  formalism  in  this  case  includes:  (i)  dividing  the  problem  into 

two  parts  -  the  (structural)  screened  structure  constant  part  and  the  potential  function  part24,25 

-  provides  a  natural  way  to  use  "frozen"  potentials'^  and  Harris-Foulks  functional^^ 
approximations  in  order  to  avoid  the  problem  of  subtracting  the  huge  total  energy  values  to 
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find  small  defect  energies.  We  note  that  frozen  potentials  have  been  used  successfully  for 

calculating  the  energetics  of  intrinsic  faults  in  Ag  using  an  Ising  model^^).  (ii)  The  influence 
of  environmental  changes  (chemical  or  structural)  on  the  structural  part  can  also  be  taken  into 

account  by  real  space  calculations28.37^ 

Unlike  previous  supercell  band  stracture^®*^^  and  layer  KKR  approaches^^>^^,  the  R- 
TBLMTO  method  may  provide  us  with  the  ability  to  simulate  (i)  impurities  in  stacking  faults 
by  going  beyond  the  single  site  approximation  (as  done  in  [38])  and,  (ii)  the  influence  of  a 
finite  concentration  of  impurities,  or  a  cluster  of  impurities  (for  treating  the  effects  of  local 
fluctuations  of  the  concentration). 

We  believe  that  developing  real  space  methods  for  studying  aluminide  materials  is  one 
possible  way  to  consider  faults  more  complex  than  planar  faults,  such  as  dislocations  which 
are  characterized  by  a  pronounced  reduction  of  crystal  symmetry  which  is  incompatible  with 
traditional  supercell  methods.  The  precision  and  teal  capabilities  of  the  proposed  R-TBLMTO 
method  were  verified  for  the  case  of  simple  planar  faults  determined  previously  by  other 
methods. 

If  aU  goes  well,  we  hope  to  use  this  technique  to  study  the  influence  of  impurities 
(i.e.,  the  chemical  effects  of  alloying)  on  the  changes  of  of  fault  energies  in  TiAl  -  including 
relaxations  at  the  faults.  This  method  appears  to  allow  us  to  calculate  self-consistently  several 

hundred  atoms  in  a  cluster  containing  several  thousand  atoms^S  and  so  may  well  provide  a 
basis  for  the  theoretical  description  of  the  faults  and  the  influence  of  impurities  on  those 
energy  parameters. 
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